Neutrino oscillations in matter and in twisting magnetic fields 
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We find the solution to the Dirac equation for a massive neutrino with a magnetic moment 
propagating in background matter and interacting with the twisting magnetic field. In frames of 
the relativistic quantum mechanics approach to the description of neutrino evolution we use the 
obtained solution to derive neutrino wave functions satisfying the given initial condition. We apply 
the results to the analysis of neutrino spin oscillations in matter under the influence of the twisting 
magnetic field. Then on the basis of the yielded results we describe spin-flavor oscillations of Dirac 
neutrinos that mix and have non-vanishing matrix of magnetic moments. We again formulate the 
initial condition problem, derive neutrinos wave functions and calculate the transition probabilities 
for different magnetic moments matrices. The consistency of the obtained results with the quantum 
mechanical treatment of spin-flavor oscillations is discussed. We also consider several applications 
to astrophysical and cosmological neutrinos. 
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I. INTRODUCTION 

Neutrino conversions from one flavor to another com- 
bined with the change of the particle helicity, e.g. ^ 
ly^, are usually called neutrino spin-flavor oscillations 
(see Ref. This neutrino oscillations scenario is im- 

portant since it could be the explanation of the time vari- 
ability of the solar neutrino flux (see, e.g., Ref. [^). Mas- 
sive flavor neutrinos are known to mix and can have non- 
zero magnetic moments. The influence of the strong mag- 
netic fleld with the realistic profile could lead to the spin- 
flavor oscillations of solar neutrinos (see, e.g., Ref. [3,]). 
Moreover, studying neutrino spin-flavor oscillations hap- 
pening inside the Sun, one will be able to discriminate 
between different solar models Q • However it was found 
out in Ref. Q that neutrino spin-flavor oscillations in so- 
lar magnetic fields give a sub-dominant contribution to 
the total conversion of solar neutrinos. 

In this paper we study neutrino spin and spin-flavor 
oscillations in matter and in an external magnetic field. 
We suppose that a neutrino is a Dirac particle with a 
non-zero magnetic moment. It should be mentioned that 
in spite of the recent claims of the experimental con- 
firmation that neutrinos are Majorana particles @, the 
question about the nature of neutrinos is still open Q. 
The possibility to distinguish between Dirac and Majo- 
rana particles in the partially polarized solar neutrino 
flux, due to the spin-flavor precession, was examined in 
Ref. i. 

To describe the evolution of the neutrino system we 
apply the technique based on the relativistic quantum 
mechanics. We start from the exact solution to the Dirac 
equation in an external field and then derive the neutrino 
wave functions satisfying the given initial condition. We 
used this method to describe neutrino flavor oscillations 
in vacuum Q , in background matter and spin-flavor 
oscillations in an external magnetic field [TTj. Note that 



neutrino spin-flavor oscillations in electromagnetic fields 
of various configurations were examined in Refs. [l^ . [isl . 
14] using the standard quantum mechanical approach. 

In Sec. HI] we find the solution to the Dirac equation 
for a neutrino propagating in background matter and 
interacting with the twisting magnetic field. Then we 
formulate the initial condition problem and receive the 
transition probability for spin oscillations in the given 
external fields. The standard quantum mechanical tran- 
sition probability formula is re-derived and the condi- 
tions of its validity are analyzed. In Sec. lIIII we apply the 
obtained Dirac equation solutions to the description of 
neutrino spin-flavor oscillations in the twisting magnetic 
fleld. First we discuss magnetic moment matrices of neu- 
trinos in flavor and mass eigenstates bases. Then we 
solve the initial condition problem in two different cases 
of the magnetic moments matrix in the mass eigenstates 
basis with (i) great diagonal elements and (ii) great non- 
diagonal elements. Note that the analogous magnetic 
moments matrices were discussed in Ref. We get 

neutrinos wave functions and calculate transition prob- 
abilities for processes like v^. The consistency 
of the Dirac-Pauli equation approach with the standard 
quantum mechanical treatment of spin-flavor oscillations, 
based on the Schrodinger evolution equation, is consid- 
ered in Sec. IIVI Then in Sec. |V] we present some appli- 
cations and finally we summarize our results in Sec. IVII 



II. NEUTRINO SPIN OSCILLATIONS IN 
MATTER AND IN A TWISTING MAGNETIC 
FIELD 

In this section we obtain the exact solution to the 
Dirac-Pauli equation for a neutrino interacting with 
background matter and a twisting magnetic fleld and dis- 
cuss spin oscillations of a single Dirac neutrino in the 
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given external fields. 

A neutrino is taken to have the non-zero mass m and 
the magnetic moment /z. The Lagrangian for this system 
has the form, 

£ = D{ij^d^ - m),y - D^luf^ - l^Da^^vF^^^ (2.1) 

where 7^ = 7^(1 + -i^)/2, a^^ = (i/2)(7^7^ - 7^7^) and 
F^^ = (E, B) is the electromagnetic field tensor. In the 
following we will discuss the situation when only mag- 
netic field B is presented, i.e. E = 0. The neutrino 
interaction with matter is characterized by the four vec- 
tor For the non- moving and unpolarized matter one 
can take that the spatial components of the vector 
are zero, i.e. f = 0. If, for instance, we consider an elec- 
tron neutrino propagating in matter, which consists of 
electrons, protons and neutrons, we obtain for the time 
component, of the vector (see, e.g., Ref. [l5l|). 

/=e,p,n 

qf={ll{^-2Q^f'^s\u'ew + 5ef), (2.2) 

where is the number density of background particles, 
J3L is the third isospin component of the matter fermion 
/, Q^f'' is its electric charge, 9w is the Weinberg angle 
and Gf is the Fermi constant. 

It should be noted that Eqs. (|2.ip and (|2.2p consti- 
tute the phenomenological model studied in the present 
paper. These expressions are valid in a relatively weak 
external magnetic field. For example, one has to take 
into account the spatial components of the vector 
if we describe neutrino propagation in background mat- 
ter composed of electrons under the influence of a very 
strong magnetic field with a/IbJ ^ max(me, T, |p|), 
where rUe is the electron mass, T is the temperature of 
background matter, 971 is its chemical potential and p is 
the neutrino momentum. This situation was analyzed in 
Ref. [3. 

Using Eq. (|2.ip one writes down the Dirac equation 
which accounts for the neutrino interaction with matter 
and magnetic field, 

ii> =T-Lv, 

H ={ap) + f3m- m/?(SB) + f{l + j^)/2, (2.3) 

where a = 7^7, /3 = 7" and S = j'^j^j are the Dirac 
matrices. Let us discuss the case of the twisting magnetic 
field, B = i?(0, sinwx, coscjx), where lo is the frequency 
of the magnetic field rotation. Sometimes it is called 
the spiral undulator magnetic field. Note that neutrino 
oscillations in twisting magnetic fields in frames of the 
quantum mechanical approach were studied in Ref. [3 | . 

We notice that the Hamiltonian Ti, in Eq. (|2.3p depends 
on neither y nor z coordinates. Therefore we assume 
that the wave function depends on these coordinates ex- 
ponentially, V ~ expiipyU + ipzz), where Py and Pz are 



constant values. Then for simplicity one can take that 
Py = Pz — 0. It means that a neutrino moves along the 
undulator axis. Let us express the neutrino wave func- 
tion in terms of the two component spinors, ly"^ = {(p, x). 
On the basis of Eq. (|2.3p we receive equations for the two 
component spinors, 

i^ =(m - ^Iia^B) + f/2)ip + {a.p, - f/2)x, 

iX =(-TO + ti{(T^B) + f/2)x + (aip, - /72)^, (2.4) 

where Px = —idx and ct±_ = {02, 0^3). 

Now we replace the neutrino wave function v with the 
new one, u ^ v — U'^v, where U = diag(il, il) and 
il — cos{iox/2) + ifji s\n.{ujx/2). Then we again express 
the new wave function using the two component spinors, 
— rj), with tp = and x = ^V- With help of the 
following properties of the matrix il: U^(cr^B)il — a^B, 
dU/dx = icriwil/2 and ilWiil — ai, as well as using 
Eq. (|2.4p we arrive to the equations for the new two com- 
ponent spinors, 

ii={m-fiBa3 + f/2)^ 

+ [{u - f )/2 + aWxH 
il) ={-m + ^iBa^ + f /2)r] 

+ [iu-f)/2 + a,px]^. (2.5) 

We notice that Eq. (|2.5p do not contain the dependence 
on X coordinate. Thus one gets that the new wave func- 
tion depends on x as ~ exp(ipa;), where p is a constant 
value, the analog of the particle momentum. It means 
that we can replace Px p in Eq. (|2.5p . 

We look for stationary solutions to Eq. (|2.5p . i.e. 
i> ~ exp(—iEt) . Supposing that this equation has a non- 
trivial solution we receive the energy levels in the form, 
E = /0/2 ± The function depends on the 

momentum and the characteristics on the external fields 
as 

E'-<^ = ^/M^ + m^+p^ - 2Ci?2, (2.6) 

where B? = \Jp'^M.'^ ■\- {^B)'^m'^ and M = 
y/{p,B)^ + /o)V4. In Eq. IIH) C = ±1 is the 
discrete quantum number. 

Using energy spectrum (|2.6p we can reproduce the re- 
sults of the previous works where the Dirac equation for 
a neutrino interacting with various external fields was 
solved. Namely, 

• neutrino interaction with a constant transversal 
magnetic field (see, e.g., Ref. [HI). This situation 
corresponds to w = and /° = 0. Using Eq. (|2.6p 

we get E = ± +p'^ — Cp-B^ that coincides 

with the energy spectrum used in Ref. [ll|; 

• neutrino interaction with background matter (see 
Ref. fl?']). This case corresponds to a; = and 
B = 0. With help of Eq. ((2?6)) we receive that E = 
f /2±yJ{p- CPJ2Y + m? that coincides with the 
results of Ref. [l3|. 
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Note that, if we set w = and i? ^ in Eq. (|2.6p . we ar- 
rive to the case of a neutrino propagating in background 
matter under the influence of a constant transversal mag- 
netic field. 

The basis spinors m^^-* and w^''^ corresponding to the 
signs ± in the dispersion relation can be found from 
Eq. (j2.5p . The general expressions for these spinors, 
which account for the particle mass exactly, are rather 
complicated. Therefore we present here the basis spinors 
for a relativistic neutrino with [m/ E) <ti 1, 



,(C) 



2^2M{M + A) 



/ fiB + CM + A' 
fiB-CM-A 
iiB-CM- A 

\fiB + CM+A 



/M ~A^CfiB\ 
yiO 1 I M-A + CfiB\ 

2J2M(M-A) A-M-CfiB 
\A-M + Ct^B I 



(2.7) 



where A = (cj — /°)/2. Note that the basis spinors in 
Eq. (|2.7p satisfy the orthonormality conditions. 



,(C)ty(C') ^ „(C)t„(C') = s^^,^ y(C)t^,(C') = 0. (2.. 



A. Neutrino evolution in matter under the 
influence of a twisting magnetic field 

Using the approach developed in our previous works [1, 
[lol . [ll| we can formulate the initial condition problem for 
the system in question. For the given initial wave func- 
tion vix^ 0) one should find the wave function ^{x, t) at 
subsequent moments of time, while a particle propagates 
in the external fields. This wave function has the form 
(see Refs. ia,iilll)), 



where 



dp 
2^' 



■'^^S{p,mp,Q), (2.9) 



+ 00 



v(p,Q)^ I dxe-'P''U\x)v{x,G), 



(2.10) 



is the Fourier transform of the initial condition for the 
fermion v and 



C=±i 



;y(C) ^y(C)t) exp (-iS('^)t) 



exp (+ 



(2.11) 



is the analog for the Pauli-Jourdan function for a spinor 
field interacting with matter and a twisting magnetic 
field. The basis spinors m'^''^ and u^''-' are presented in 
Eq. (1131) • To derive Eqs. (l^ - ipUl) we use orthonor- 
mality of the basis spinors 



Let us suppose that initially a neutrino is in the state 
with the following wave function: v{x, 0) = e''^^^0; where 
= (1/2)(1, -1,-1,1). It is possible to check that 
(1/2)(1 — Si)^o — <^o- Hence, the spinor u[x,Q) de- 
scribes a particle propagating along the a;-axis, with its 
spin directed opposite to the a;-axis, i.e. a left-handed 
neutrino. Analogous initial condition was adopted in 
Refs. [§, [l^, [ll| where neutrino flavor and spin-flavor 
oscillations were studied. 

Using Eq. fTTOl) we find that v{p, 0) = 2tt5{p - k - 
ll}/2)^q. It is interesting to note that the following iden- 



tity is satisfied: (■ 



,(0 



® V 



(C)tl 



^0 = 0. Therefore no par- 



ticles with "negative" energies appear in neutrino inter- 
acting with considered external fields. Using Eqs. (|2.7p 
and (|2.9p - (|2.1ip as well as the chosen initial condition we 
arrive to the right-polarized component of the final wave 
function, 



i/^(x,t) =i(l + Si)Ka;,i) 

= exp[i(fc + w)2:-i/°i/2] 



(2.12) 



^lB 

2M 



where = (1/2)(1, 1, 1, 1). 

Supposing that initially no right-polarized particles are 
present and with help of Eq. (|2.12[) we calculate the tran- 
sition probability for the process ^ 



" (/iB)2 + Z\2 

2 fE+-E- 



X sm 



(2.13) 



p=fc+w/2 



It can be seen from Eq. (|2.13|1 that the resonance in neu- 
trino spin oscillations occurs when Z\ — > 0. One finds 
from Eq. (HH) that {E+ - £'")/2 = -fiB at zi = 0. 
Therefore the resonance transition probability is always 
Prcs(i) =sm^ (nBt). 

To analyze Eq. (|2.13p we introduce the group velocity. 



V(C) 



dE 



dp EiO 



1-C 



i?2 



(2.14) 



Now we can distinguish three different cases. 



1. First we suppose that p ~ 0. This situation can 
happen if fc ~u;/2. With help of Eq. 

we obtain that the energy levels are 

E(0 = 

\J(jn — C,[iBY + A?- Using Eq. (|2.14p we receive 
that the group velocity vanishes, V^''-' = 0. It means 
that the neutrino is captured by the twisting mag- 
netic field. 

2. Now we assume that p = ±A^, with p 7^ 0. For the 
dcfiniteness we discuss the situation when p = M. 
since the case p = —M. can be considered anal- 
ogously. The energies corresponding to different 
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values of ( are 



E- =2M 1 



(2.15) 



Using Eq. (|2.14p one can compute the group veloc- 
ities, 



V 



2M- 



4_21-l/2 



V" =1 



1 + 3 



7W2 2M^ 
{p.Bf^ 



(2.16) 



In Eqs. (|2.15|1 and (|2.16p we suppose that M. 
On the basis of Eqs. ([^1^ and ([^1^ we get the 
resonance energies [A 0), 



V2hB' 
1 - 



res 

and group velocities 



A{^,By 



^ V- ^1- 

V2' 2(//i3)2 



(2.17) 



(2.18) 



It should be noted that group velocities are always 
less than one, [see, e.g., Eq. (|2.18|) ]. 

3. The last situation is realized when p ^ ±A^ and 
p 7^ 0. The energies in this case have the form. 



2ip-CM) 



1-C 



ipBl 
Mp 



(2.19) 



The expression for the transition probability (|2.13p 
is now rewritten in the following way: 



Pit) 



if^Bf 



ifiB)^ + Z\2 

X sin^ (y^lBpT^t) 



(2.20) 



Note that transition probability expressions for 
spin oscillations derived earlier (see, e.g., Ref. [T3 |) 
coincide with Eq. (|2.20p which is valid only if 
p ±M and p ^ 0. 

It should be noted that the "non-standard" regimes 
in neutrino spin oscillations described in items [T] and [2] 
are likely to be realized for neutrinos with small initial 
momenta (see also Sec. |V] below). 



III. NEUTRINO SPIN-FLAVOR OSCILLATIONS 
IN A TWISTING MAGNETIC FIELD 

Now we apply the results of the previous section to 
the description of neutrino spin-flavor oscillations in a 
twisting magnetic field. Let us study the evolution of 
two Dirac neutrinos (i^a, vp) that mix and interact with 
the external electromagnetic field F^^. The Lagrangian 
for this system has the form 



\=a,p X\'=a,p 



mx\'V\V\' 



(3.1) 



Here (toaa') and {M\\i) are the mass and the mag- 
netic moments matrices that are generally independent. 
By definition these matrices are intoduced in the flavor 
eigenstates basis. The electromagnetic field is taken to 
have the same configuration as in Sec. [Til 

To analyze the dynamics of the system we again set the 
initial condition by specifying the initial wave functions 
of the flavor neutrinos v\ and then analytically determine 
the field distributions at following moments of time. We 
assume that the initial condition is 



i'q(x,0) = 0, vp{x,G) ^ S,{x), 



(3.2) 



where ^(x) is a function to be specified. One of the pos- 
sible choices for the initial condition for vp is the plane 



wave field distribution, £,{x) 



'ikx 



Co (see Refs. (9|, M, 



[Ti]|). If we study ultrarelativistic initial particles, we can 
choose the spinor as in Sec. [TTl i.e. in the following 
form: = (i/2)(i, _i, _l, l). 

In order to eliminate the vacuum mixing term in 
Eq. p.ip . i.e. to diagonalize the mass matrix, we intro- 
duce a new basis of the wave functions, the mass eigen- 
state basis ipa, a = 1,2, obtained from the original flavor 
basis i>\ through the unitary transformation 



(3.3) 



where the matrix {U\a) is parametrized in terms of a 
mixing angle as usual 



COS f — sm ( 



I sin 6 COS t 



(3.4) 



The Lagrangian (|3.ip rewritten in terms of the fields 
takes the form 



-^(■01, ■02) = ^ 'Co(V'a) 
a=l,2 

^ ^ Pab'ipaCTtJ.ulpbF'^'' , (3.5) 



afc=l,2 
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where £,o{ipa) = '^ai^l^d^ — ma)'4'a is the Lagrangian for 
the free fermion i/jq with the mass nia and 



(3.6) 



\\'=a,l3 



is the magnetic moment matrix presented in the mass 
eigenstates basis. Using Eqs. (I3.2I) - (|3.4I) the initial con- 
ditions for the fermions Tpa become 



i^i{x, 0) ~ sm9^{x), i'2{x, 0) = cos0^(x) 



(3.7) 



For the given configuration of the electric and magnetic 
fields we write down the Dirac-Pauli equation for ^l)a^ 
resulting from Eq. (13. 5p . as follows: 



(3.8) 



\'4)a^'Hai}a^Vi}b, a, 5=1,2, a^b, 



where Tia — (ap) + (3raa — /iQ/3(SB) is the Hamilto- 
nian for the particle accounting for the magnetic field, 
V = — /i/3(SB) describes the interaction of the transi- 
tion magnetic moment with the external magnetic field, 
A*a = Maa, arid ^ — /ii2 — /i2i are elements of the matrix 

(/iab)- 

To find the general solution to Eq. (|3.8p we follow the 
method used in Sec.HTland introduce the new wave func- 
tions i/'a — ^^^a- All the calculations are identical to 
those made in Sec. [Til Therefore we present the final 
result for the wave functions V'a, 



V'a(a;,0 



+00 



dp 

7^ 



X ^ [ai^)(i)^i^)exp(-i£;(';)0 
C=±i 

+ 5(0(t)^;(C)exp(+ii?(0i)\ 
.(0 



(3.9) 



where the energy levels E)t' are 
Here [see Eq. (EH)] 



Ma ^^{tiaBY+uJ-'/^. 



(3.10) 



(3.11) 



The basis spinors up and vP can be obtained from 



Eq. (|2.7p by the following replacement: n ^ iia, -M. ^ 
Aia and /° — > 0. Our main goal is to determine the 
non-operator coefficients aP and bP so that to satisfy 



both the initial condition (|3.7[) and the evolution equa- 
tion (|3.8|) . Generally the coefficients aP (t) and bP (t) 
arc functions of time. 



Spin-flavor oscillations in case of diagonal 
magnetic moments 



fia ^ A*- This case should be analyzed with help of the 
perturbation theory. We expand the wave functions "00 
in a series 



^a(x,t)=^(°) (a:, 0+^-^(0^,*) + 



(3.12) 



where ^^°'>{x,t) corresponds to the solution of Eq. 
when we neglect the potential V there. The function 
ipP{x,t) is linear in the transition magnetic moment /i 
etc. We omit terms of higher order in fi in Eq. (|3.12p . 
They can be also accounted for but the corresponding 
calculations arrear to be cumbersome in the general case. 

Using orthonormality conditions of the basis spinors 
[see also Eq. 



and the results of our previous work [ll| (see also Sec. HI]) 
we can receive from Eq. (13.91) the expression for the zero 

order (in ^) wave functions ipP , which correspond to the 
first term in Eq. ([XT^ . 



^P{x,t)^U{x) 
where 



%^P-Saip,t)MP,0)- (3.13) 

ZTT 



+ 00 



dxe-'P''U''{x)'tl;a{x,0), (3.14) 



is the Fourier transform of the initial condition for the 
spinor V'a- Here 



Saip,t)= J2 

C=±i 



+ («P®«F)exp(+ii?Pi) 



(3.15) 



is the analog of the Pauli-Joudan function in the twisting 
magnetic field [see also Eq. (|2.1ip ]. 

Using Eqs. (|3.13p - (|3.15p for the given initial condition 
we can find the wave functions at any subsequent mo- 
ments of time. For example, if one initially has the left- 
handed neutrino i^p , then field distribution of the right- 
handed component of the fermion i>a is 



1 



^) = 2 (1 + ^1 ) [cos Oi;i {x,t)- sin 6*^2 {x, t)] 



: sm t) cos ue 



i(A:+a; )x 



(3.16) 



X 



2M 

2M2 



-iE7t 



Kq. 



p=k+uj/2 



In this section we suppose that magnetic moments ma- 
trix in the mass eigenstates basis is close to diagonal, i.e. 



To receive Eq. p.l6p we use the same technique as in 
Sec. ini Therefore we may omit the details of calcula- 
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tions. On the basis of Eq. (|3.16p one obtains the transi- 
tion probabihty for the process vj^ — > in the form 



P 



(0) 



At) 



sin ^l^t 



Ml 



— — sm ( 1 

Ml 



4——— sm 
M1M2 



E+ - 



X sm I — —t 



X sm 



El + E^ - E2 - E2 



(3.17) 



Note that the phase of oscillations in Eq. (|3.20p de- 
pends on the frequency of the twisting magnetic field. 
It should be noted that, if we put a; = in Eqs. p.l9p 
and p.20p , the transition probability coincides with that 
from our work [ll[ where we studied neutrino spin-flavor 
oscillations in the constant transversal magnetic field. 

If we studied the special case of massive neutrinos hav- 
ing equal magnetic moments, /ii = = ^J■o, we would 
obtain the expected result from Eq. (|3.19p . Namely, 
Eq. (|3.19p can be rewritten as P = PpPs, where Pp ~ 
sin^(20) sin^[<I>(fc)t] is the usual transition probability of 
fiavor oscillations and 



Ps = —7^2 — sm [ilst), 



(3.21) 



The energies in Eqs. (|3.16p and (|3.17p are given 

in Eq. (|3.10p . In Eq. (|3.17p we also suppose that p — 
k + Lo/2. 

The analysis of Eq. (|3.17p is almost identical to that in 
Secini Therefore we present in the explicit form the final 
results for the wave function and the transition probabil- 
ity in the most important case when p max(ma, Ma)- 
This situation corresponds to spin-flavor oscillations of 
ultrarelativistic neutrinos. Now the wave function of i>a 
becomes 



zy^^-*^ (x, t) ==i sin cos exp[i(fc -I- uj)x — ipt] 

riiB f .m'i \ . . , 
—— exp -i-r-t smMit 
Mi\2p) 

^exp(^-r^.jsmM,t 



(3.18) 



p=k+u/2 



and the transition probability in Eq. (j3.17D is 



(0) _sin^ 



HiB . fi2B . 

sm Alii — — sinM2t 



Ml 
M1M2 



(3.19) 

2 



M2 

sinTWit sin7W2i sin^[$(fc)] L 



In Eq. p.l9p we use the notation for the oscillations 
phase. 



$(fc) 



dm? 



A{k + ul2) 



(3.20) 



and Sni? = — ra\ is the mass squared difference. In 
deriving Eqs. p.lSp and p.l9p we use the analog of the 
energy expansion in Eq. (j2.19p . 



is the probability of spin oscillations between differ- 
ent polarization states within each mass eigenstate. In 
Eq. (IX^ ns = ^/{jM0fTW/2^. That is, since the 
magnetic moment interactions are insensitive to fiavor, 
the transitions between flavors are solely due to the mass 
mixing. 

One can obtain the first order corrections (linear in /i) 
to Eqs. (|3.18p and (|3.19p . These corrections correspond 
to the second term in Eq. p.l2p . The expressions for the 
corrections to the mass eigenstates wave functions are 



V^i'Hx,t) 



27r 



(3.22) 



iU{:x) 

J —00 

C=±i 

(v'^P (g> v^p^) exp {+iEi'^h)vn^p] Mp, 0), 



where 



CjK) = f dt'exp{+iEi<h')S,{p,t'), 
Jo 



dt' exp (-iEi<h') Ship, t') 



(3.23) 



In Eqs. ((3?22)) and ((3?23l) a 7^ 6. For the details of the 
derivation of Eqs. (|3.22p and p.23p the reader is referred 
to Ref. 11] and Sec. [Til of the present paper. Note that 
V = -fiBpj:3 in Eq. ([2221) • 

The calculations of the the first order corrections based 
on Eqs. p.22p and p.23p are rather cumbersome. There- 
fore we present here only the final results in the case when 
p 3> max(ma, A^a). One has the expression for the cor- 
rection to the wave function. 
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i{k-\-uj)x 



1 



4:MlM2 

- {fiiH2B^ - M1M2 ~ t^V4) cos 20 
+ {M1-M2)' 



{HiH2B^ +M1M2 - cos 20 



sin^i sinAt 



d 



A 



In Eq. (|3.24p we use the notations, 



D 



{M1+M2)- 



sinrft _ sin-P^ i5f 



(3.24) 



E+ + E+ 



El 



E = 



2 

E^ + 



p + T(/c) -^A^, 

:^p+T(fc)+5M, 



and 



A 



E+ - E^ 



E^ - E2 



^{k)-SM, d = 
Q>{k) + 6M, D 



E+ - E- 



E^ — E2 



$(fc) 
a $(fc) 



where 



Tik) 



m\ + 7712 

4(fc + w/2)' 



5M = , M = . 



To obtain Eq. p.24p we use the identity {v^\v\£,q) = 0, which means that no antineutrinos are produced. 

On the basis of Eqs. p.lSp and (|3.24p one calculates the correction to the transition probability which has the 
following form: 



(1) _ nBsin2e 



(^71^2^2 +Mi7W2 -c^V4) 



cos 26* 

Zi 



$(fc) sin[2$(fc)t 
6M 



f Ub B Us B 
' siiiMitcosM2t — — cos Alii sinAl2i 



''^^sin^Xit + ^sin^X^t- 



Mi 



M2 



M2 

^J■lB ^ H2B 

Ml M2 



sinTMii sinA^2i 



2 + sinAlitsinAl2isin2[<i>(fc)t]^ 



{tiUi2B^ - MiM2-oj^/^) 



cos 26* 
Z2 



M 
2 



uiB . , . ^ 
' sin^Alit 



Ml 

( _ A^2^ 

\Mi M2 



^sin^M2t 
M2 



$(fc) sin[2$(fc)t] 
/ y.iB 112B 



uiB U2B 

smMitcosM2t r-^ cos A^it sinA^2t 



V Ml M2 



Ml 

sinA^it sin A^2i 



M2 



+ iMi-M2) 



$(fc) 



2Zi 
Mi-B • 2 



sin Alit sin Al2isin^ [^ik)t] 



6Msin\2^(k)t] { sinMitcosM2t - cos Mit sin M2t 
* Ml M2 



mB . 



M 



fiiB M2-B\ . 



- sin^ A^ii + ain^ M2t - ^-r— + ^-r— sin A^ii sin A^2i 
1 M2 \Mi M2 ) 



2 + 4^ ^ sin7Witsin7W2tsin2[$(fc)f 

\Mi M2 ) 



-{M1+M2) — 



7Wsin[2$(fc)i] 



Uj B Li B 

— — sii\MitcosM2t r— cos Alii sinAl2* 

All At2 
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^k)( ^ sin' M^t-!^ sin' M2t+ ( - f^) sinM^tsinM^t 
\Mi M2 \Mi M2J 

2 - sin Mit sin M2tsin'mk)t]^ |, 



(3.25) 



where 



h\^:\t) in Eq. dSl 



(3.26) 



Note that, if we again put a; = in Eqs. (|3.24p and (|3.25l) . 
we reproduce the results of our previous work [11]. 

It can be noticed from Eqs. (|3.25p and p.26p that the 
perturbative approach is vahd until Z12 7^ 0. If cither 
Z\ or Z2 is equal to zero, we can expect that some non- 
perturbative effects like resonances can occur. Unfortu- 
nately these effects cannot be quantitatively described in 
frames of the approach based on the Dirac-Pauli equa- 
tion used in the present work. To analyze such phenom- 
ena one should carry out numerical computations within 
the Schrodinger equation approach (see also Sec. lIVp . 
Nevertheless we evaluate the possibility that Zi = Q for 
spin-flavor oscillations between active and sterile neutri- 
nos (see, e.g., Ref. in the twisting magnetic field 
of the Sun. In this case one can take into account the 
magnetic moment of an active neutrino only. For the fol- 
lowing parameters: k lOMeV, Sm' ^ 10"^ eV^ [11], 
A*activ ~ 10-"jUB, B - 100 kG and w - IQ-i^eV [H, we 
can see that the quantities /iactiv^, <i>(fc) and to are of the 
same order of magnitude of 10~^^ eV. Thus the violation 
of the validity of the perturbation theory is quite possible 
for this kind of situation and resonance phenomena can 
happen. 

The sum of Eqs. p.l9p and p.25p gives one the tran- 
sition probability of spin-flavor oscillations up to terms 
linear in fj, in case of the magnetic moments matrix which 
is close to diagonal. 



B. Spin-flavor oscillations in case of non-diagonal 
magnetic moments 



In this section we study neutrino spin-flavor oscilla- 
tions in case of the non-diagonal magnetic moments ma- 
trix. It means that the transition magnetic moment 
dominates over the diagonal ones, i.e. we assume that 
^J■ > ^J■a■ 

Now we start directly from Eq. (|3.9p . However one can- 
not treat the potential V as the small perturbation. To 
solve this problem we should use the method elaborated 
in Ref. [ill] . The following ordinary differential equations 
can be derived to determine the coefficients al^\t) and 



Y: [apu(<\.p{-iEPt) 



ib^P 



C'=±i 
: exp (— 



exp {-iE^'^'k) 



C'=±i 

+ bpv^<">exp{+iEl'^'h) 



(3.27) 



where it+ — kq, = = (1/2)(1, 1, —1, —1) and 

v^"^ = (1/2)(1, — 1, 1, — 1). The quantum number ( is 
the eigenvalue of the operator (Sp)/]p] = Ei, i.e. for 
example, Eiw^''^ = Cw^^'- Note that the current definition 
of C differs from that we use in previous sections. We 
also drop the subscripts a and b in the basis spinors since 
we study the evolution of ultrarelaticistic neutrinos and 
assume that <C H- The energies Ei^'^ in Eq. p.27p 
take the form 



Ei^ = ^ml + {p + Cio/2r. (3.28) 

For the details of the derivation of Eq. (|3.28p and the 
basis spinors the reader is referred to Sec. [iTl of this work. 

The initial conditions should be added to the differen- 
tial equations (|3.27p . 

4«(o) = -^^.(^)t^.(p,o), 



(3.29) 



2tt 



Eq. p.29p results from Eq. (|3.9p and the orthonormality 
of the basis spinors m*^^) and v^''\ 

Taking into account the following identi- 
ties: (m^^^JF ]«(«')) = 0, {u^\V\u^) = and 
(M^jyju^) = —I^B, which can be verified by means of 
direct calculations, one reveals that Eq. (|3.27p is reduced 
to the form 



i = -af^iBexp [iiE^-Ej^)t] 



(3.30) 

(C) 



Note that the analogous equation for the functions ba 
can be also obtained from Eq. (|3.27p . Eq. (|3.30l) is similar 
to that analyzed in Ref. [llj] (see also Ref. pj)]). There- 
fore we write down its solution, e.g., for the functions 



"^2(0 = i^sinfi±texp(±iw±t/2)aj;2(0), (3.31) 
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where 



n± = y/ilJ-B)"^ + (w±/2)2, uj± = 2^{k)±uj. (3.32) 



In deriving Eq. (j3.3ip we take into account that initiaUy 
only left-handed neutrinos are presented, i.e. a^(0) = 0. 
Indeed V'a(p, 0) ~ and with help of Eq. p.29p we get 

thata+(0) = 0. 

Finally, using Eqs. ([23]), (03), (EH) and ([?3T|) we ar- 
rive to the right-handed component of I'a , 



(3.33) 



X exp 



1 P + 



4p 



9 „sin Q+t , 
X I cos 6 — sm^ ( 



sin fl_t 



Kq. 



With help of Eq. (|3.33p we can compute the transition 
probability for the process like vj^ —> in case of mag- 
netic moments matrix with great non-diagonal elements, 

(3.34) 



n „siiiQ+t 9„sinr2_t 
cos'^ - snr 0—— 



1 2 



Note that, if we approach to the limit w = in Eq. (|3.34p . 
we reproduce the result of our work where spin- 
flavor oscillations of neutrinos with similar magnetic mo- 
ments matrix were studied. 

Let us analyze neutrino oscillations at small fre- 
quencies of the twisting magnetic field, co — 
\J inBy + [<I>(A;)]2. In this situation the transition prob- 
ability in Eq. ()3.34|) can be rewritten as 



P{t) =A{t)sm^{not), 

A{t) ^A^,n, + 2SAsm^{5nt). 



(3.35) 



where SA — {Am, 
w$(fc)/(2r2o) and 



Amin)/2, 6n = {n+ - o_)/2 



An 



COS 29 



2 r 



1 -V^ COS 26* 



"o 



COS 26* 



(3.36) 



Eqs. p.35p and (|3.36p show that the behavior of the sys- 
tem is analogous to beatings occurring in interference of 
two oscillations with different amplitudes and frequen- 
cies. 

Let us discuss neutrinos with the following parameters: 
sin^e = 0.3, 6m^ = IQ-^eV^ [11,^= IQ-^Vb and 
k = 100 MeV. It is known that rather strong twisting 
magnetic fields, up to the critical value of i? ~ 10^'' G, 
can exist in the early Universe [20| . The time depen- 
dence of neutrino oscillations probability is schematically 
depicted on Fig. [T] for such a magnetic field strength and 
u! = 10""'^'^ eV. As one can see on this figure, the rapidly 




t,s 



FIG. 1; The time dependence of the neutrino oscillations 
probability in the twisting magnetic field in the case fi ^ fii^2 
sit small values of lu. 



varying transition probability P{t) (solid line) is modu- 
lated by the slowly varying function A{t) (dashed line). 
This time dependence is different from that described in 
Ref. [13. 

It is also possible to see on Fig. [T] that the typical 
time scale of the amplitude modulation of the transition 
probability is about Ta ~ 0.1s. The production rate of 
right-handed neutrinos in the early Universe should be 
less than the expansion rate of the Universe in order not 
to affect the primordial nucleosynthesis [2l[. Hence one 
has Ta/i > 1, where 



1.24 X 10 



3 I 



V 100 MeV 



is the Hubble parameter [22| and Tpi is the primordial 
plasma temperature. Supposing that neutrinos are at 
thermal equilibrium at 100 MeV, i.e. k ~ Tp\, we get 
that TAh - 100 > 1. 

Despite the mentioned above discrepancy it is inter- 
esting to compare the result of this section [Eq. (|3.34p ] 
with the analogous transition probability formula for 
Majorana neutrinos at small mixing angle 9^0. 
Note that in this situation magnetic moments matrices 
in Eqs. ([33]) and (031) [o r Eq. jSJ])] coincide. In this 
limit we obtain from Eq. p.34p the following transition 
probability: 



P,l^,r(0 =^^sin2(0Mt), 



(3.37) 



where = v^(^i3)2 + [<i>(fc) + cj/2]^. It can be seen 
that Eq. (|3.37p coincides with the transition probabilities 
derived in Ref. [14] where spin-fiavor oscillations of Ma- 
jorana neutrinos in twisting magnetic fields were studied 
on the basis of the quantum mechanical approach. 
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It should be also noticed that the transition proba- 
bility in Eq. p.34p vanishes at high frequencies, w » 
ma.x{nia, fiB), due to the dependence of the oscillations 
phase on oj [see Eq. (13.201) 1. This phenomenon was also 
mentioned in our paper 13] in which we examined spin- 
flavor oscillations of Majorana neutrinos in rapidly vary- 
ing external fields. 



Then we look for the solutions to the Schrodinger 
equation idvE'/dt = i/^E', with the Hamiltonian given in 
Eq. (|4.5p . in the form ^ ^ e^''^*. The secular equation 
for A is the forth order algebraic equation in general case. 
However it can be solved in two situations. 



A. Diagonal magnetic moments matrix 



IV. QUANTUM MECHANICAL DESCRIPTION 
OF NEUTRINO OSCILLATIONS IN A 
TWISTING MAGNETIC FIELD 

In this section we demonstrate that the analog of the 
main results, yielded in Sees. IIIlXl and ITlI Bl within the 
Dirac-Pauli equation approach, can be also obtained with 
help of the standard quantum mechanical treatment of 
spin-flavor oscillations. The consistency between these 
two approaches is discussed. 

To study neutrinos evolution in frames of the 
Schrodinger equation approach it is convenient to make 
the coordinates transformation. We assume that 
k = {0,0, k) and B = B{smujt,cosu!t,0). Now the 
Schrodinger equation and the effective Hamiltonian for 
the neutrinos mass eigenstates have the form, 



i* = H = 



H 



mass 

t 

B 



Hb 



(4.1) 



where iJmass = diag(£i,f2), ^1,2 



+ k'^ and 



Hb — — {Bx + iBy 



= -\{nab)Be~ 



(4.2) 



where /ii.2 and are the elements of the magnetic mo- 
ments matrix (|3.6p . The neutrinos wave function has the 
following form: "^"^ = (V']", j "0?; V'?)? where V'1'2^ are 
one-component functions. It can be seen that Eqs. (|4.ip 
and (|4.2p is the generalization, for the case of Dirac neu- 
trinos, of the corresponding expressions used in Ref. p^ . 
The initial condition for the wave function follows from 
Eq. (j3Jl) . 



^-^(0) = (sin 0, cos 6*, 0,0). 
Let us make the matrix transformation, 



m =diag(e-'"*/^ e~'^'/^, e'^''^, e''^*/^) 



(4.4) 



The Hamiltonian H governing the time evolution of the 
modified wave function 4' is presented in the form 



H 



i{Hab)B Hmass + tuj/2 



(4.5) 



where 1 is the 2x2 unit matrix. The initial condition 
for ^1/ coincides with that for 4" [see Eq. (|4.3p ] due to the 
special form of the matrix QJ in Eq. 



In the case when ^ = and ^1,2 7^ the roots of the 
secular equation are 



A 



8a±Ma, a =1,2, 



(4.6) 



where A^a is defined in Eq. ()3.1ip . The basis spinors 
, which are the eigenvectors of the Hamiltonian H are 
expresed in the following way: 



,+ - 




(4.7) 



where Z1.2 — \J M.\^2 + w/2. Note that the vectors 
correspond to the eigenvalues A^ . 

The general solution to the Schrodinger evolution 
equation has the form, 

E E 4'^^^i'^)exp(-iAi';)i), (4.8) 

a=l,2C=±l 

(4.3) where the coefficients a'S' should be chosen so that to 
satisfy the initial condition in Eq. (|4.3p . We choose these 
quantities as 

1 i^i-S _ . /) 2i 



: sm ( 



, 1 iM25 



V2M2 Z2 



V2Mi 



(4.9) 



Then using Eqs. (l44l) and (l46 l) - ((49l) we get the right- 
polarized components of the wave function ^ in the form 



ipf{t) =exp[-i(£i -uj/2)t] sm0sm{Mit) 
Tpf{t) =exp[-i(£2 - w/2)i] cos6'sin(A^2t) 



fiiB 
Ml ' 

M2^ 

Mo ■ 



(4.10) 
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Finally taking into account Eqs. (|3?3)) . ((O)) and (|4?T0l) 

wc arrive to the right-handed component of Va, 



form 



Juit/2 



exp(-if2^) sin{M2t) 



exp(— i£ii) sin(A^it) 



M2 



Ml 
(4.11) 



One can see from Eq. (14. lip that the expression for 
obtained in frames of the Schrodinger approach coincides 
(to within some irrelevant phase factor) with the analo- 
gous expression derived using the Dirac-Pauli equation 
[see Eq. ([XT^ ]. 



B. Non-diagonal magnetic moments matrix 

In the situation when /ii.2 = and /i 7^ the secular 
equation can be also solved analytically and the corre- 
sponding roots are 



(4.12) 



ipf{t) ^nB cos 9 exp[-i(£ - Lu/2)t] ^"^(^"+^) 



-I- 

sin(ri_i) 



■02 (i) =Ai^sin6'exp[-i(£: - uj/2)t] — ^ 



(4.16) 



On the basis of Eq. (|4.16p we receive the wave function 



yR (t) ^^iB cxp[-i{£ -uj/2)t] 



, „sin Q+t . , 
X I cos 9 — —— sm^ 



sin f2_i 



(4.17) 



Comparing Eq. (j4.17p with the analogous expres- 
sion (|3.33p derived in frames of the Dirac-Pauli equation 
approach wc again find an agreement to within the phase 
factor. 

It should be however mentioned that using the quan- 
tum mechanical treatment of spin-flavor oscillations one 
cannot reproduce the expression for the phase of the neu- 
trino oscillations (|3.20p . It means that in Eqs. (|4.1ip 
and (|4.17p we have the standard quantum mechanical 
vacuum oscillations phase $QA/(fc) = Sm? /{Ak). 



where £ — {£1 + £2)72 and Vl± are given in Eq. (|3.32p . 
The eigenvectors of the Hamiltonian H have the following 
form: 



■"1,2 



Via = 




(4.13) 



where R± = \/^+ ± ^+/2, S± — yj^- ± uj-/2 and ijJ± 
are given in Eq. (I3.32p . The spinors ui^ and vi^ in 
Eq. (|4.13p correspond to the eigenvalues X^^ and 
respectively. 

The general solution to the Schrodinger equation for 
the function ^I^ takes the form, 

^{t) — exp(— i£i)[ai'Ui exp(— ifi+t) + q:2'U2 exp(ir2+t) 
-I- Pivi exp(-ir2_t) + (32V2 exp(if7_i)], (4.14) 

We again choose the coefficients ai^ and (ii^ in 
Eq. (|4.14p to satisfy the initial condition in Eq. (|4.3p . 
These coefficients have to be chosen as 

ai2=±icos6' — ^ /9i2=±isin6' — , (4.15) 



With help of Eqs. (g^) and ((TT^ - (|i?T5)) we obtain the 
right-handed components of the wave function ^' in the 



V. APPLICATIONS 

Let us discuss the applicability of our results to one 

specific oscillation channel, vj;^ ^ v^. According to the 
recent neutrino oscillations data (see, e.g., Ref. [19|) the 
mixing angle between and Vr is close to its maximal 
value of 7r/4. For such a mixing angle the magnetic mo- 
ment matrix given in Eq. p.6p is expressed in the form. 



[Mr 



+ M^^]/2 + Mr 



-[Mrr-M^^]/2 



- \Mrr - M^,^\ 12 [Mrr + M^^]/2 - Mr, 



(5.1) 



Eq. p.lOp is valid when this matrix is close to diago- 
nal, i.e. if \{Mrr - M^^)/2\ < \{Mrr + M^^)/2 ± Mr^,\. 
In contrast to the mixing angles, experimental data and 
theoretical predictions for the values of neutrino mag- 
netic moments are not very reliable [23j . However it is 
known that the diagonal magnetic moments M\x could 
be very small in the extensions of the standard model 
Mxx ~ 10~^^(?7iAA/eV)/iB [111- The transition magnetic 
moments, M^fi in our notations, can be much greater up 
to the experimental limit of 10^^"/iB [25]. One can see 
that for any conceivable values of the masses of the known 
neutrinos, -/V/^^ and M^t are several orders of magnitude 
smaller than 10^^°/iB. Our result (|3.19p is valid in this 
case. 

It is worth mentioning that Eqs. p.l3p and (I3.22p are 
in principle applicable for particles with arbitrary initial 
conditions, e.g., with small initial momenta. Therefore 
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one can discuss the evolution and oscillations of relic neu- 
trinos. These particles can gravitationally cluster in the 
Galaxy [2^. It is possible to relate the temperatures of 
relic neutrinos and CMB protons, and respectively, 

= (^-j T, « 0.72T^. (5.2) 

For the present value = 2.7° K we get k. 
1.93° K [26|. Using the estimate for the neutrino mass 
m ~ 0.1 eV, because the sum of all neutrinos masses 
should be less than 1 eV [27[ , and taking into account 
that these neutrinos are non-relativistic particles, one 
obtains the typical momentum of a relic neutrino k ^ 
7 X 10"'^ eV. For example, to realize the "non-standard" 
neutrino propagation regime described in item[T]of SecUll 
one should use an undulator with the frequency |a;| ~ 2fc 
or with the period L = 27r/w ~ 0.1 mm. 

Strong periodic electromagnetic fields, with a short 
spatial oscillations length, can be found in crystals [28j . 
In Ref. [2PI it was proposed to manufacture deformed 
crystals with submillimeter periods. The undulator radi- 
ation was recently reported in Ref. ^ to be produced in 
undulators with periods of 0.1—1 mm. Therefore artificial 
crystalline undulators with required periods are used in 
various experiments and hence they can serve as a good 
tool to explore properties of rehc neutrinos. It should be 
mentioned that neutrino scattering on a polarized tar- 
get and possible tests of neutrino magnetic moment were 
examined in Ref. Isil. 



VI. SUMMARY 

We have described the evolution of Dirac neutrinos in 
matter and in a twisting magnetic field. We have applied 
the recently developed approach (see Refs. d, [lO, ITlj ) 
which is based on the exact solutions to the Dirac equa- 
tion in an external field with the given initial condition. 

First (Sec. |TT| we have found the solution to the Dirac 
equation for a neutral 1/2-spin particle weakly interact- 
ing with background matter, that is equivalent to an ex- 
ternal axial-vector field, and non-minimally coupled to 
an external electromagnetic field due to the possible pres- 
ence of an anomalous magnetic moment. We have dis- 
cussed the situation when a neutrino interacts with the 
twisting magnetic field. The energy spectrum and ba- 
sis spinors have been obtained. We have applied these 
results to derive the transition probability of spin os- 
cillations in matter under the influence of the twisting 
magnetic field. The scope of the standard quantum me- 
chanical approach to the description of neutrino spin os- 
cillations has been analyzed. 

Then (Sec. IIIip we have used the obtained solution to 
the Dirac equation for the description of neutrino spin- 
flavor oscillations in a twisting magnetic field. We sup- 
posed that two Dirac neutrinos could mix and have non- 
vanishing matrix of magnetic moments. Moreover the 



mass and magnetic moments matrices in the flavor eigen- 
states basis are generally independent, i.e. the diagonal- 
ization of the mass matrix, that means the transition to 
the mass eigenstates basis, does not lead to the diago- 
nalization of the magnetic moments matrix. We have 
discussed two possibilities. 

In Sec, nil Al we have assumed that magnetic moments 
matrix in the mass eigenstates basis has great diagonal 
elements compared to the non-diagonal ones. In this case 
one can analyze neutrino spin-flavor oscillations pertur- 
batively. Note that the perturbative approach allows one 
to discuss neutrinos with an arbitrary initial condition. 
For instance, the evolution of particles with small initial 
momenta can be accounted for. The appearance of non- 
perturbative phenomena like resonances is analyzed with 
an example of active-to-sterile neutrinos oscillations. We 
have discussed the opposite situation in Sec. lIII I3l i.e. the 
magnetic moment matrix with the great non-diagonal el- 
ements. In this case one had to treat the evolution of the 
system non-perturbatively. We have demonstrated that 
this situation is analogous to beatings resulting from the 
superposition of two oscillations. In both cases we have 
obtained neutrino wave functions, consistent with the ini- 
tial conditions, and the transition probabilities. Note 
that all the results are in agreement with our previous 
work (ll| if we set uj — Q, i.e. discuss a constant transver- 
sal magnetic field. We have also examined some limiting 
cases and compared our results with the previous studies. 

It has been shown in Sec. IIVI that one can derive 
the analog of the major results obtained in Sees. IIII Al 
and IIII Bl using the Schrodinger evolution equation ap- 
proach for the description of spin-flavor oscillations of 
Dirac neutrinos. The correspondence between these two 
approaches has been considered. 

In Sec. |V| we have discussed magnetic moments ma- 
trices in various theoretical models which predict neutri- 
nos magnetic moments. The validity of our approach for 
these situations has been considered. The applications of 
our results to the studying of cosmological neutrinos in 
laboratory conditions have been examined. In particular 
we have suggested that artificial crystalline undulators 
could be useful for such a research. 

The results obtained in the present work are valid for 
arbitrary magnetic field strength. The general case of 
spin-flavor oscillations of Dirac neutrinos in a twisting 
magnetic field with an arbitrary magnetic moment ma- 
trix has not been studied analytically earlier. Both ex- 
perimental and theoretical information about magnetic 
moments of Dirac neutrinos is known to be very lim- 
ited (see, e.g. Refs. [H, |2^). Therefore our results can 
be helpful since they enable one to describe phenomeno- 
logically spin-flavor oscillations of Dirac neutrinos under 
the influence of the magnetic field in question provided 
neutrinos possess non-vanishing matrix of magnetic mo- 
ments. Although we consider neutrinos, it is possible to 
straightforwardly apply our formalism to the description 
of any 1/2-spin particles. 



13 



Acknowledgments 

The work has been supported by the Academy of Fin- 
land under the contract No. 108875. The author is thank- 
ful to the Russian Science Support Foundation for a grant 



as well as to Efrain Ferrer (Western Illinois University) 
and Kimmo Kainulainen (University of Jyvaskyla) for 
useful comments. The referees' remarks are also appre- 
ciated. 



[1] J. Schechter and J. W. F. Valle, Phys. Rev. D 24, 1883 
(1981); M. B. Voloshin, M. I. Vysotskii, and L. B. Okun', 
JETP 64, 446 (1986); Sov. J. Nucl. Phys. 44, 440 (1986); 
E. Akhmedov, Phys. Lett. B 213, 64 (1988); C.-S. Lim 
and W. J. Marciano, Phys. Rev. D 37, 1368 (1988). 

[2] B. C. Chauhan, J. Pulido and R. S. Raghavan, JHEP 
0507, 051 (2005), hep-ph/0504069; J. Barranco et ai, 
Phys. Rev. D 66, 093009 (2002), hep-ph/0207326. 

[3] A. B. Balantekin and P. Loreti, Phys. Rev. D 48, 5496 
(1993); A. Friedland and A. Gruzinov, Astropart. Phys. 
19, 575 (2003), hep-ph/0202095. 

[4] M. Picarieno et al, to be pubUshed in JHEP, 0705.4070 
[hep-ph]. 

[5] E. Kh. Akhmedov and J. PuUdo, Phys. Lett. B 553, 7 
(2003), hep-ph/0209192; A. B. Balantekin and C. Volpe, 
Phys. Rev. D 72, 033008 (2005), hep-ph/0411148. 

[6] H. V. Klapdor-Klcingrothaus, et al., Phys. Lett. B 586, 
198 (2004), hcp-ph/0404088. 

[7] S. R. Elliott and J. Engel, J. Phys. G 30, R183 (2004), 
hep-ph/0405078. 

[8] V. B. Semikoz, Nucl. Phys. B 498, 39 (1997), hep- 
ph/9611383. 

[9] M. Dvornikov, Phys. Lett. B 610, 262 (2005), hcp- 
ph/0411101; in Proceedings of the IPM school and con- 
ference on Lepton and Hadron Physics, Tehran, 2006, ed. 
by Y. Farzan, eConf C0605151 (2007), hep-ph/0609139; 
hep-ph/0610047. 

[10] M. Dvornikov, Eur. Phys. J. C 47, 437 (2006), hep- 
ph/0601156; to be published in J. Phys. Conf. Ser., 
0708.2975 [hep-ph]. 

[11] M. Dvornikov and J. Maalampi, Phys. Lett. B 657, 217 
(2007), hep-ph/0701209; M. Dvornikov, to be published 
in proceedings of the 14*'' International Baksan School 
"Particles and Cosmology", 0708.3572 [hep-ph]. 

[12] M. S. Dvornikov and A. I. Studonikin, Phys. At. Nucl. 
64, 1624 (2001); tbtd. 67, 719 (2004); E. Kh. Akhmedov 
and M. Yu. Khlopov, Sov. J. Nucl. Phys. 47, 689 (1988); 
A. Egorov, A. Lobanov, and A. Studenikin, Phys. Lett. 
B 491, 137 (2000), hep-ph/9910476. 

[13] M. Dvornikov, Phys. At. Nucl. 70, 342 (2007), hep- 
ph/0410152. 



[14] A. Yu. Smirnov, Phys. Lett. B 260, 161 (1991); E. Kh. 
Akhmedov, P. I. Krastev, and A. Yu. Smirnov, Z. Phys. 
G 48, 701 (1991); E. Kh. Akhmedov, S. T. Petcov, and 
A. Yu. Smirnov, Phys. Rev. D 48, 2167 (1993), hep- 
ph/9301211. 

[15] M. Dvornikov and A. Studenikin, JHEP 09, 016 (2002), 

hep-ph/0202113. 
[16] E. EUzalde, E. J. Ferrer and V. de la Incera, Phys. Rev. 

D 70, 043012 (2004), hep-ph/0404234. 
[17] A. Studonikin and A. Ternov, Phys. Lett. B 608, 107 

(2005), hep-ph/0412408; A. E. Lobanov, Phys. Lett. B 

608, 136 (2005), hep-ph/0506007. 
[18] P. Keranen et al., Phys. Lett. B 597, 374 (2004), hep- 

ph/0401082. 

[19] M. Maltoni et al, New J. Phys. 6, 122 (2004), hep- 
ph/0405172. 

[20] H. Athar, Phys. Lett. B 366, 229 (1996). 

[21] A. D. Dolgov, Phys. Rcpt. 370, 333 (2002), hep- 

ph/0202122; M. Dvornikov, Int. J. Mod. Phys. A 21, 

2403 (2006), hep-ph/0501063. 
[22] S. Weinberg, Gravitation and Cosmology: Principles and 

Applications of the General Theory of Relativity (Wiley, 

New York, 1972), p. 534. 
[23] N. F. Bell et al, Phys. Rev. Lett. 95, 151802 (2005), 

hep-ph/0504134; N. F. Bell et al., Phys. Lett. B 642, 

377 (2006), hep-ph/0606248. 
[24] K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 45, 963 

(1980); M. Dvornikov and A. Studenikin, Phys. Rev. D 

69, 073001 (2004), hcp-ph/0305206. 
[25] W.-M. Yao et al., J. Phys. G 33, 1 (2006). 
[26] A. Ringwald and Y. Wong, JGAP 0506, 013 (2005), hep- 

ph/0410136. 

[27] S. Hannestad and G. G. Raffelt, JGAP 0611, 016 (2006), 

astro-ph/0607101. 
[28] U. I. Uggerh0j, Rev. Mod. Phys. 77, 1131 (2005). 
[29] S. Bcllucci et al, Phys. Rev. Lett. 90, 034801 (2003), 

physics/0208028. 
[30] V. T. Baranov et al, JETP Lett. 82, 562 (2005). 
[31] T. I. Rashba and V. B. Semikoz, Phys. Lett. B 479, 218 

(2000), hep-ph/0003099. 



